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INTRODUCTION 

Aquatic  invertebrates  are  aptly  applied  to  bioassessment  since  they  are  known  to 
be  important  indicators  of  stream  ecosystem  health  (Hynes  1970).  Long  lives,  complex 
life  cycles  and  limited  mobility  mean  that  there  is  ample  time  for  the  benthic  community 
to  respond  to  cumulative  effects  of  environmental  perturbations. 

This  report  summarizes  data  collected  in  August  2001  from  3  sites  on  the 
Stillwater  River,  Park  County,  Montana.  Aquatic  invertebrate  assemblages  were  sampled 
by  personnel  of  the  Montana  Department  of  Environmental  Quality  (DEQ).  Study  sites 
lie  within  the  Middle  Rockies  ecoregion.  (Woods  et  al.  1999).  A  multimetric  approach  to 
bioassessment  such  as  the  one  applied  in  this  study  uses  attributes  of  the  assemblage  in  an 
integrated  way  to  measure  biotic  health.  A  stream  with  good  biotic  health  is  "  . . .a 
balanced,  integrated,  adaptive  system  having  the  full  range  of  elements  and  processes  that 
are  expected  in  the  region's  natural  environment. . ."  (Karr  and  Chu  1999).  The  approach 
designed  by  Plafkin  et  al.  (1989)  and  adapted  for  use  in  the  State  of  Montana  has  been 
defined  as  "...  an  array  of  measures  or  metrics  that  individually  provide  information  on 
diverse  biological  attributes,  and  when  integrated,  provide  an  overall  indication  of 
biological  condition."  (Barbour  et  al.  1995).  Community  attributes  that  can  contribute 
meaningfully  to  interpretation  of  benthic  data  include  assemblage  structure,  sensitivity  of 
community  members  to  stress  or  pollution,  and  functional  traits.  Each  metric  component 
contributes  an  independent  measure  of  the  biotic  integrity  of  a  stream  site;  combining  the 
components  into  a  total  score  reduces  variance  and  increases  precision  of  the  assessment 
(Fore  et  al.  1995).  Effectiveness  of  the  integrated  metrics  depends  on  the  applicability  of 
the  underlying  model,  which  rests  on  a  foundation  of  three  essential  elements  (Bollman 
1998).  The  first  of  these  is  an  appropriate  stratification  or  classification  of  stream  sites, 
typically,  by  ecoregion.  Second,  metrics  must  be  selected  based  upon  their  ability  to 
accurately  express  biological  condition.  Third,  an  adequate  assessment  of  habitat 
conditions  at  each  site  to  be  studied  is  needed  to  assist  in  the  interpretation  of  metric 
outcomes. 

Implicit  in  the  multimetric  method  and  its  associated  habitat  assessment  is  an 
assumption  of  correlative  relationships  between  habitat  parameters  and  the  biotic  metrics, 
in  the  absence  of  water  quality  impairment.  These  relationships  may  vary  regionally, 
requiring  an  examination  of  habitat  assessment  elements  and  biotic  metrics  and  a  test  of 
the  presumed  relationship  between  them.  Bollman  (1998)  has  recently  studied  the 
assemblages  of  the  Montana  Valleys  and  Foothill  Prairies  ecoregion,  and  has 
recommended  a  battery  of  metrics  applicable  to  the  montane  ecoregions  of  western 
Montana.  This  metric  battery  has  been  shown  to  be  sensitive  to  impairment,  related  to 
habitat  assessment  parameters,  and  consistent  over  replicated  samples. 

Habitat  assessment  enhances  the  interpretation  of  biological  data  (Barbour  and 
Stribling  1991),  because  there  is  generally  a  direct  response  of  the  biological  community 
to  habitat  degradation  in  the  absence  of  water  quality  impairment.  If  biotic  health  appears 
more  damaged  than  the  habitat  quality  would  predict,  water  pollution  by  metals,  other 
toxicants,  high  water  temperatures,  or  high  levels  of  organic  and/or  nutrient  pollution 
might  be  suspected.  On  the  other  hand,  an  "artificial"  elevation  of  biotic  condition  in  the 
presence  of  habitat  degradation  may  be  due  to  the  paradoxical  effect  of  mild  nutrient  or 
organic  enrichment  in  an  oligotrophic  setting. 


METHODS 

Aquatic  invertebrates  were  sampled  by  Montana  DEQ  personnel  on  August  2  and 
24  at  3  sites  on  the  Stillwater  River.  Site  locations  and  sampling  dates  are  indicated  in 
Table  1 .  The  sampling  method  employed  was  that  recommended  in  the  Montana 
Department  of  Environmental  Quality  (DEQ)  Standard  Operating  Procedures  for  Aquatic 
Macroinvertebrate  Sampling  (Bukantis  1998).  In  addition  to  aquatic  invertebrate  sample 
collection,  habitat  quality  was  visually  evaluated  at  each  site  and  reported  by  means  of 
the  habitat  assessment  protocols  recommended  by  Bukantis  (1998). 

Evaluated  habitat  features  include  instream  conditions,  larger-scale  channel 
conditions  including  flow  status,  streambank  condition,  and  extent  of  the  riparian  zone. 
Scores  were  assigned  in  the  field  to  each  habitat  measure,  and  these  scores  were  totaled 
and  compared  to  the  maximum  possible  score  to  give  an  overall  assessment  of  habitat. 

Aquatic  invertebrate  samples  and  associated  habitat  data  were  delivered  to 
Rhithron  Biological  Associates,  Missoula,  Montana,  for  laboratory  and  data  analyses. 
In  the  laboratory,  the  Montana  DEQ-recommended  sorting  method  was  used  to  obtain 
subsamples  of  at  least  300  organisms  fi"om  each  sample,  when  possible.  Organisms  were 
identified  to  the  lowest  possible  taxonomic  levels  consistent  with  Montana  DEQ 
protocols. 

To  assess  aquatic  invertebrate  communities  in  this  study,  a  multimetric  index 
developed  in  previous  work  for  streams  of  western  Montana  ecoregions  (Bollman  1998) 
was  used.  Multimetric  indices  result  in  a  single  numeric  score,  which  integrates  the 
values  of  several  individual  indicators  of  biologic  heaUh.  Each  metric  used  in  this  index 
was  tested  for  its  response  or  sensitivity  to  varying  degrees  of  human  influence. 
Correlations  have  been  demonstrated  between  the  metrics  and  various  symptoms  of 
human-caused  impairment  as  expressed  in  water  quality  parameters  or  instream, 
streambank  and  stream  reach  morphologic  features.  Metrics  were  screened  to  minimize 
variability  over  natural  environmental  gradients,  such  as  site  elevation  or  sampling 
season,  which  might  confound  interpretation  of  results  (Bollman  1998).  The  multimetric 
index  used  in  this  report  incorporates  multiple  attributes  of  the  sampled  assemblage  into 
an  integrated  score  that  accurately  describes  the  benthic  community  of  each  site  in  terms 
of  its  biologic  integrity.  In  addition  to  the  metrics  comprising  the  index,  other  metrics, 
which  have  been  shown  to  be  applicable  to  biomonitoring  in  other  regions  (Kleindl  1995, 
Patterson  1996,  Rossano  1995)  were  used  for  descriptive  interpretation  of  results.  These 
metrics  include  the  number  of  "dinger"  taxa,  long-lived  taxa  richness,  the  percent  of 
predatory  organisms,  and  others.  They  are  not  included  in  the  integrated  bioassessment 
score,  however,  since  their  performance  in  western  Montana  ecoregions  is  unknown. 
However,  the  relationship  of  these  metrics  to  habitat  conditions  is  intuitive  and 
reasonable. 

The  six  metrics  comprising  the  bioassessment  index  used  in  this  study  were 
selected  because  both  individually  and  as  an  integrated  metric  battery,  they  are  robust  at 
distinguishing  impaired  sites  fi-om  relatively  unimpaired  sites  (Bollman  1998).  In 
addition,  they  are  relevant  to  the  kinds  of  impacts  that  are  present  in  the  Stillwater  River 
drainage.  They  have  been  demonstrated  to  be  more  variable  with  anthropogenic 
disturbance  than  with  natural  environmental  gradients  (Bollman  1998).  Each  of  the  six 
metrics  developed  and  tested  for  western  Montana  ecoregions  is  described  below. 


Table  1.  Sampling  sites  and  dates:  three  sites  on  the  Stillwater  River.  Sites  are  listed  from 
upstream  to  downstream. 

^'^^                     Site  description                   Sampling                             GPS 
designation Date Latitude Longitude 

SROl  above  Daisy  Creek  8/2/01  45°  03' 59"    109°  59' 38" 

SR02  below  Daisy  Creek  8/2/01  45°  04*  06"     109°  59' 39" 

SR04  near  wilderness  boundary  8/24/01  45°  05' 33"     109°  59*  32" 


1.  Ephemeroptera  (mayfly)  taxa  richness.  The  number  of  mayfly  taxa 
declines  as  water  quality  diininishes.  Impairments  to  water  quality  which  have  been 
demonstrated  to  adversely  affect  the  ability  of  mayflies  to  flourish  include  elevated 
water  temperatures,  heavy  metal  contamination,  increased  turbidity,  low  or  high  pH, 
elevated  specific  conductance  and  toxic  chemicals.  Few  mayfly  species  are  able  to 
tolerate  certain  disturbances  to  instream  habitat,  such  as  excessive  sediment 
deposition. 

2.  Plecoptera  (stonefly)  taxa  richness.  Stoneflies  are  particularly 
susceptible  to  impairments  that  affect  a  stream  on  a  reach-level  scale,  such  as  loss  of 
riparian  canopy,  streambank  instability,  channelization,  and  alteration  of 
morphological  features  such  as  pool  frequency  and  function,  riffle  development  and 
sinuosity.  Just  as  all  benthic  organisms,  they  are  also  susceptible  to  smaller  scale 
habitat  loss,  such  as  by  sediment  deposition,  loss  of  interstitial  spaces  between 
substrate  particles,  or  unstable  substrate. 

3.  Trichoptera  (caddisfly)  taxa  richness.  Caddisfly  taxa  richness  has  been 
shown  to  decline  when  sediment  deposition  affects  their  habitat.  In  addition,  the 
presence  of  certain  case-building  caddisflies  can  indicate  good  retention  of  woody 
debris  and  lack  of  scouring  flow  conditions. 

4.  Number  of  sensitive  taxa.  Sensitive  taxa  are  generally  the  first  to 
disappear  as  anthropogenic  disturbances  increase.  The  list  of  sensitive  taxa  used  here 
includes  organisms  sensitive  to  a  wide  range  of  disturbances,  including  warmer  water 
temperatures,  organic  or  nutrient  pollution,  toxic  pollution,  sediment  deposition, 
substrate  instability  and  others.  Unimpaired  streams  of  western  Montana  typically 
support  at  least  four  sensitive  taxa  (BoUman  1998). 

5.  Percent  filter  feeders.  Filter-feeding  organisms  are  a  diverse  group;  they 
capture  small  particles  of  organic  matter,  or  organically  enriched  sediment  material, 
from  the  water  column  by  means  of  a  variety  of  adaptations,  such  as  silken  nets  or 
hairy  appendages.  In  forested  montane  streams,  filterers  are  expected  to  occur  in 
insignificant  numbers.  Their  abundance  increases  when  canopy  cover  is  lost  and 
when  water  temperatures  increase  and  the  accompanying  growth  of  filamentous  algae 
occurs.  Some  filtering  organisms,  specifically  the  Arctopsychid  caddisflies 
(Arctopsyche  spp.  and  Parapsyche  sp.)  build  silken  nets  with  large  mesh  sizes  that 
capture  small  organisms  such  as  chironomids  and  early-instar  mayflies.  Here  they  are 


considered  predators,  and,  in  this  study,  their  abundance  does  not  contribute  to  the 
percent  filter  feeders  metric. 

6.    Percent  tolerant  taxa.  Tolerant  taxa  are  ubiquitous  in  stream  sites,  but 
when  disturbance  increases,  their  abundance  increases  proportionately.  The  list  of 
taxa  used  here  includes  organisms  tolerant  of  a  vsdde  range  of  disturbances,  including 
warmer  water  temperatures,  organic  or  nutrient  pollution,  toxic  pollution,  sediment 
deposition,  substrate  instability  and  others. 

Scoring  criteria  for  each  of  the  six  metrics  are  presented  in  Table  2.  Metrics  differ 
in  their  possible  value  ranges  as  well  as  in  the  direction  the  values  move  as  biological 
conditions  change.  For  example,  Ephemeroptera  richness  values  may  range  fi-om  zero  to 
ten  taxa  or  higher.  Larger  values  generally  indicate  favorable  biotic  conditions.  On  the 
other  hand,  the  percent  filterers  metric  may  range  from  0%  to  100%;  in  this  case,  larger 
values  are  negative  indicators  of  biotic  health.  To  facilitate  scoring,  therefore,  metric 
values  were  transformed  into  a  single  scale.  The  range  of  each  metric  has  been  divided 
into  four  parts  and  assigned  a  point  score  between  zero  and  three.  A  score  of  three 
indicates  a  metric  value  similar  to  one  characteristic  of  a  non- impaired  condition.  A  score 
of  zero  indicates  strong  deviation  fi-om  non-impaired  condition  and  suggests  severe 
degradation  of  biotic  health.  Scores  for  each  metric  were  summed  to  give  an  overall 
score,  the  total  bioassessment  score,  for  each  site  in  each  sampling  event.  These  scores 
were  expressed  as  the  percent  of  the  maximum  possible  score,  which  is  18  for  this  metric 
battery. 

Table  2.  Metrics  and  scoring  criteria  for  bioassessment  of  streams  of  western  Montana 
ecoregions  (Bollman  1998). 


Score 

metric 

3 

2 

1 

0 

Ephemeroptera  taxa  richness 

>5 

5-4 

3-2 

<2 

Plecoptera  taxa  richness 

>3 

3-2 

1 

0 

Trichoptera  taxa  richness 

>4 

4-3 

2 

<2 

Sensitive  taxa  richness 

>3 

3-2 

1 

0 

Percent  filterers 

0-5 

5.01  -  10 

10.01-25 

>25 

Percent  tolerant  taxa 

0-5 

5.01-10 

10.01-35 

>35 

The  total  bioassessment  score  for  each  site  was  expressed  in  terms  of  use-support. 
Criteria  for  use-support  designations  were  developed  by  Montana  DEQ  and  are  presented 
in  Table  3a.  Scores  were  also  translated  into  impairment  classifications  according  to 
criteria  outlined  in  Table  3b. 

In  this  report,  certain  other  metrics  were  used  as  descriptors  of  the  benthic 
community  response  to  habitat  or  water  quality  but  were  not  incorporated  into  the 
bioassessment  metric  battery,  either  because  they  have  not  yet  been  tested  for  reliability 
in  streams  of  western  Montana,  or  because  results  of  such  testing  did  not  show  them  to  be 
robust  at  distinguishing  impairment,  or  because  they  did  not  meet  other  requirements  for 
inclusion  in  the  metric  battery.  These  metrics  and  their  use  in  predicting  the  causes  of 
impairment  or  in  describing  its  eflFects  on  the  biotic  community  are  described  below. 


The  modified  biotic  index.  This  metric  is  an  adaptation  of  the  Hilsenhoflf  Biotic 
Index  (HBI,  HilsenhofF  1987),  which  was  originally  designed  to  indicate  organic 
enrichment  of  waters.  Values  of  this  metric  are  lowest  in  least  impacted 
conditions.  Taxa  tolerant  to  saprobic  conditions  are  also  generally  tolerant  of 
warm  water,  fine  sediment  and  heavy  filamentous  algae  growth  (Bollman, 
unpublished  data).  Loss  of  canopy  cover  is  often  a  contributor  to  higher  biotic 
index  values.  The  taxa  values  used  in  this  report  are  modified  to  reflect  habitat 
and  water  quality  conditions  in  Montana  (Bukantis  1998).  Ordination  studies  of 
the  benthic  fauna  of  Montana's  foothill  prairie  streams  showed  that  there  is  a 
correlation  between  modified  biotic  index  values  and  water  temperature,  substrate 
embeddedness,  and  fine  sediment  (Bollman  1998).  In  a  study  of  reference 
streams,  the  average  value  of  the  modified  biotic  index  in  least-impaired  streams 
of  western  Montana  was  2.5  (Wisseman  1992). 

Taxa  richness.  This  metric  is  a  simple  count  of  the  number  of  unique  taxa  present 
in  a  sample.  Average  taxa  richness  in  samples  fi-om  reference  streams  in  western 
Montana  was  28  (Wisseman  1992).  Taxa  richness  is  an  expression  of  biodiversity, 
and  generally  decreases  with  degraded  habitat  or  diminished  water  quality. 
However,  taxa  richness  may  show  a  paradoxical  increase  when  mild  nutrient 
enrichment  occurs  in  previously  oligotrophic  waters,  so  this  metric  must  be 
interpreted  with  caution. 

Percent  predators.  Aquatic  invertebrate  predators  depend  on  a  reliable  source  of 
invertebrate  prey,  and  their  abundance  provides  a  measure  of  the  trophic 
complexity  supported  by  a  site.  Less  disturbed  sites  have  more  plentifiil  habitat 
niches  to  support  diverse  prey  species,  which  in  turn  support  abundant  predator 
species. 

Number  of  "dinger"  taxa.  So-called  "dinger"  taxa  have  physical  adaptations  that 
allow  them  to  cling  to  smooth  substrates  in  rapidly  flowing  water.  Aquatic 
invertebrate  "dingers"  are  sensitive  to  fine  sediments  that  fill  interstices  between 
substrate  particles  and  eliminate  habitat  complexity.  Animals  that  occupy  the 
hyporheic  zones  are  included  in  this  group  of  taxa.  Expected  "dinger"  taxa 
richness  in  unimpaired  streams  of  western  Montana  is  at  least  14  (Bollman, 
unpublished  data). 

Number  of  long-lived  taxa.  Long-lived  or  semivoltine  taxa  require  more  than  a 
year  to  completely  develop,  and  their  numbers  decline  when  habitat  and/or  water 
quality  conditions  are  unstable.  They  may  completely  disappear  if  channels  are 
dewatered  or  if  there  are  periodic  water  temperature  elevations  or  other 
interruptions  to  their  life  cycles.  Western  Montana  streams  with  stable  habitat 
conditions  are  expected  to  support  six  or  more  long-lived  taxa  (Bollman, 
unpublished  data). 


Table  3a.  Criteria  for  the  assignment  of  use-support  classifications  /  standards  violation 
thresholds  (Bukantis,  1997). 

%  Comparability  to  reference 

Use  support 

>75 

25-75 

<25 

Full  support—standards  not  violated 

Partial  support—moderate  impairment- 
standards  violated 

Non-support-severe  impairment— standards 
violated 

Table  3b.  Criteria  for  the  assignment  of  impainnent  classifications  (Plafkin  et  al.  1989). 

%  Comparability  to  reference 

Classification 

>83 

54-79 
21-50 
<17 

nonimpaired 
slightly  impaired 
moderately  impaired 
severely  impaired 

RESULTS 

Habitat  assessment 

Figure  1  compares  habitat  assessment  results  for  the  3  sites  in  this  study.  Table  4 
itemizes  the  evaluated  habitat  parameters  and  shows  the  assigned  scores  for  each. 

Figure  1.  Total  habitat  assessment  scores  for  sites  on  the  Stillwater  River,  August  2001. 
Sites  are  listed  from  upstream  to  downstream. 


OS 

s 


100 
90 
80 
70 
60 
50 
40 
30 
20 
10 
0 


^^^H^t- 

.    H^. 

^^^P^»^  >■  '* 

1 1 

' 

SROl 


SR02 


SR04 


Table  4.  Stream  and  riparian  habitat  assessment.  Stillwater  River.  August  2001.  Not  all 
parameters  were  evaluated  at  site  SR04;  total  score  is  based  only  on  those  parameters 
actually  assessed. 


Max.  possible 
toore 

Parameter 

SROl 

SR02 

SR04 

10 

Riffle  development 

10 

10 

10 

10 

Benthic  substrate 

5 

3 

7 

20 

Embeddedness 

19 

6 

15 

20 

Channel  alteration 

20 

20 

20 

20 

Sediment  deposition 

20 

12 

20 

20 

Channel  flow  status 

16 

12 

18 

20 

Bank  stability 

8/8 

8/8 

10/10 

20 

Bank  vegetation 

8/8 

6/6 

n.a. 

20 

Vegetated  zone 

2/2 

9/3 

n.a. 

160 

Total 

126 

103 

110 

Percent  of  maximum 

79 

64 

92 

CONDITION* 

OPT 

SUB 

OPT 

♦Condition  categories:  Optimal  (OPT)  >  80%  of  maximum  score;  Sub-optimal  (SUB)  75  -  56%;  Marginal 
(MARG)  49  -  29%;  Poor  <23%.  Adapted  from  Plafkin  et  al.  1988. 

Habitat  assessment  scores  indicate  optimal  habitat  conditions  at  the  upper  and 
lower  sites  (SROl  and  SR04)  and  sub-optimal  conditions  at  the  intermediate  site  (SR02). 
At  site  SROl,  benthic  substrate  diversity  was  perceived  to  be  marginal;  it  was  described 
as  dominated  by  gravel,  but  with  other  particle  sizes  present  too.  Riparian  zone  width  was 
judged  poor. 

At  the  intermediate  site  (SR02),  marginal  benthic  substrate  diversity  was  again 
noted,  although  the  presence  of  several  size  categories  of  particles  were  reported  in  field 
notes.  Embeddedness  was  perceived  to  be  extensive.  The  width  of  the  riparian  zone  was 
judged  marginal  on  one  side  of  the  channel  but  optimal  on  the  other  side. 

At  the  downstream  site  (SR04),  diversity  of  substrate  sizes  was  not  perceived  to 
be  a  problem,  and  embeddedness  was  judged  slight.  Streambank  vegetative  protection 
and  riparian  zone  width  were  not  evaluated,  since  bedrock  dominated  bank  structure. 

Bioassessment 

Macroinvertebrate  taxa  lists,  metric  results,  and  other  information  for  each  sample 
are  given  in  the  Appendix.  Figure  2  compares  the  total  bioassessment  scores  calculated 
for  macroinvertebrate  communities  collected  at  each  of  the  3  sites.  Breakdown  of  scores 
for  each  metric  calculated  fi"om  Stillwater  River  aquatic  invertebrate  samples  is  presented 
in  Table  5. 

Applying  this  method  to  the  data  generated  Irom  Stillwater  River  samples 
suggests  that  the  upper  and  intermediate  sites  (SROl  and  SR02)  fully  supported 
designated  uses  and  were  essentially  unimpaired  in  terms  of  biotic  integrity,  however,  an 
examination  of  the  taxonomic  con^sition  of  the  assemblage  collected  at  site  SR02 
suggests  severe  impairment,  possibly  associated  with  metals  contamination.  This  is 
explained  in  the  following  section.  Low  abundance  of  organisms  in  the  sample 
complicated  the  evaluation  of  the  lower  site  (SR04);  but  the  taxonomic  composition  of 


Figure  2.  Total  bioassessment  scores,  expressed  as  percent  of  maximum,  for  three  sites 
on  the  Stillwater  River.  August  2001.  The  revised  metric  battery,  used  here,  does  not 
detect  the  evidence  for  severe  metals  pollution  at  two  of  the  three  sites.  See  text  for 
explanation. 


SROl 


SR02 


SR04 


the  assemblage  suggests  that  the  presence  of  heavy  metals  influences  the  benthic 
community  at  this  site. 

Site  SROl  yielded  slightly  fewer  Plecoptera  and  Trichoptera  taxa  than  expected 
for  a  montane  stream.  At  site  SR02,  fewer  Ephemeroptera  taxa  were  collected  than 
expected.  No  Ephemeroptera  taxa  were  present  in  the  sample  taken  at  SR04. 

Aquatic  invertebrate  communities 

The  upper  Stillwater  River  site  (SROl)  supported  a  sensitive,  diverse,  cold-water 
assemblage  comprised  of  rich  mayfly,  stonefly,  and  caddis  fly  fauna.  Among  the 
sensitive  organisms  taken  in  the  collection  were  the  stoneflies  Zapada  columbiana, 
Megarcys  sp.,  and  Yoraperla  sp.;  together  these  taxa  comprised  20%  of  sampled 
organisms.  They  were  the  only  stonefly  taxa  collected,  however,  resulting  in  a  diminished 
Plecoptera  richness  score,  which  may  suggest  some  degree  of  reach-scale  habitat 
disturbance  or  inadequacy,  such  as  an  abbreviated  riparian  zone,  channel  alteration,  or 
extensive  streambank  instability.  Only  3  caddisfly  taxa  were  collected,  but  12  "dinger" 
taxa  were  present;  fine  sediment  deposition  apparently  did  not  compromise  habitats  at 
this  site.  The  low  biotic  index  score  (2.55)  and  abundant  sensitive  taxa  (5)  indicated  good 
water  quality.  The  dominant  taxon  at  the  site  was  the  midge  Diamesa  sp.,  a  cold- 
stenotherm;  it  comprised  29%  of  organisms  in  the  sample.  All  expected  functional 
components  were  amply  represented  in  the  assemblage. 

Only  4  mayfly  taxa  were  collected  at  the  intermediate  site  (SR02),  and  the 
modified  biotic  index  value  (5.38)  was  higher  than  expected  for  a  montane  site.  Impaired 
water  quality  may  be  inferred  fi-om  these  findings.  While  the  high  biotic  index  value  may 
have  been  the  result  of  a  fortuitous  collection  of  a  "bloom"  of  midges  in  the  Cricotopus 


Table  5.  Metric  values  and  bioassessments  for  the  Stillwater  River.  August  2001. 
The  assessment  classification  and  use  support  designation  in  parentheses  is  based  on  a 
sample  with  inadequate  numbers  of  organisms.  The  revised  metric  battery,  used  here, 
does  not  detect  the  evidence  for  severe  metals  pollution  at  two  of  the  three  sites.  See  text 
for  explanation. 


SITES 

SROl 

SR02 

SR04 

METRICS 

METRIC  VALUES 

Ephemeroptera  richness 

6 

4 

0 

Plecoptera  richness 

3 

5 

2 

Trichoptera  richness 

3 

6 

5 

Number  of  sensitive  taxa 

5 

9 

4 

Percent  fliterers 

1 

0 

0 

Percent  tolerant  taxa 

1 

1 

0 

METRIC  SCORES 

Ephemeroptera  richness 

3 

2 

0 

Plecoptera  richness 

2 

3 

2 

Trichoptera  richness 

2 

3 

3 

Number  of  sensitive  taxa 

3 

3 

3 

Percent  filterers 

3 

3 

3 

Percent  tolerant  taxa 

3 

3 

3 

TOTAL  SCORE  (max.=18) 

16 

17 

14 

PERCENT  OF  MAX. 

89 

94 

(78) 

Impairment  classiflcation* 

NON 

NON 

(SLD 

USE  SUPPORT  T 

FULL 

FULL 

(FULL) 

1 .  Classifications:  (NON)  non-impaired,  (SLI)  slightly  impaired,  (MOD)  moderately  impaired,  (SEV) 

severely  impaired  See  Table  3b. 

*Use  support  designations:  See  Table  3a. 


Isocladius  Group,  which  comprised  62%  of  sampled  organisms,  their  abundance,  along 
with  the  decline  in  species  richness  and  abundance  of  mayflies,  is  consistent  with 
patterns  associated  with  metals  contamination.  The  composition  of  the  benthic 
community  at  this  site,  compared  to  that  of  the  site  above  Daisy  Creek,  demonstrates  a 
shift  in  the  abundance  of  mayflies  toward  midges  in  the  genus  Cricotopus.  This  is 
illustrated  in  Figure  3.  In  addition,  heptageniid  mayflies  {Cinygmula  sp.,  Epeorus  sp.,  and 
Rhithrogena  sp),  which  accounted  for  more  than  25%  of  the  sample  collected  at  the 
upper  site  (SROl),  represented  only  3%  of  organisms  in  the  sample  collected  below  the 
confluence  with  Daisy  Creek.  The  severe  decrease  in  abundance  of  heptageniid  mayflies 
and  the  taxonomic  shifts  fi'om  abundant  mayflies  and  few  Orthocladiinae  to  few  mayflies 
and  abundant  Orthocladiinae  strongly  suggests  that  contamination  by  metals  is  present  at 
the  lower  site  (Clements  1994).  Thus  it  appears  that  the  scores  derived  fi-om  applying  the 
revised  metric  battery  do  not  adequately  assess  this  condition. 

The  lower  site  (SR04)  produced  a  sample  with  too  few  organisms  for  conclusions 
to  be  reliable,  it  is  unclear  fi-om  the  data  whether  the  sample  represents  actual  low 


abundance  of  benthos  at  the  site,  or  whether  sampling  effort  was  inadequate. 
Nonetheless,  some  observations  can  be  made  based  on  the  taxonomic  composition  of  the 
sampled  assemblage.  Mayflies  were  completely  absent,  strongly  suggesting  that  metals 
contamination  was  also  severe  here.  Among  the  dominant  taxa  was  Cricotopus  sp.  Shifts 
in  taxonomic  composition  of  benthic  communities  appear  to  persist  at  this  downstream 
she. 

CONCLUSIONS 

•  The  upper  site  (SROl)  on  the  Stillwater  River  supported  a  diverse  assemblage 
characteristic  of  an  unimpaired  montane  stream. 

•  At  the  intermediate  site  (SR02),  the  revised  metric  battery  did  not  detect  the 
strong  evidence  for  metals  contamination.  Instead,  the  likelihood  of  severe 
pollution  was  evident  in  the  taxonomic  composition  of  the  sampled  assemblage. 
Heptageniid  mayflies,  which  were  abundant  at  the  upper  site,  were  scarce  below 
the  confluence  of  Daisy  Creek.  In  addition,  the  midge  Cricotopus  sp. 
overwhelmed  the  community  there. 

•  Low  sample  abundance  makes  observations  about  the  benthic  assemblage  at  SR04 
uncertain,  but  evidence  of  metals  contamination  persist  at  this  site. 

•  Each  site  supported  only  a  single  long-lived  taxon,  which  suggests  the  possibility 
that  seasonal  loss  of  surface  flow  may  periodically  interrupt  life  cycles  in  these 
reaches  of  the  Stillwater  River. 
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Above  Daisy  Creek 


Below  Daisy  Creek 


■  Misc.  taxa 
D  Trichoptera 


Ephemeroptera  D  Plecoptera 
Diptera  ■  Chironomidae 


Figure  3.  A  comparison  of  the  relative  abundances  of  the  major  taxonomic  groups 
present  at  sites  SROl  (above  Daisy  Creek)  and  SR02  (below  Daisy  Creek),  illustrating  the 
shift  in  composition  of  assemblages  suggestive  of  severe  metals  contamination. 
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Aquatic  Invertebrate  Taxonomic  Data 


Site  Name:  Stillwater  River:  above  Daisy  Creek 


Site  ID:  Y04USR01    8/2/01 
Taxon 

Approx  pcfxxnt  of  sample  used:  40 

Quantity          Percent 

HBI 

FFG 

Acari 

5                 1.58 

5 

PA 

Total  Misc.  Taxa 

5                1.58 

Baetis  tricaudatus 
Drunella  coloradensis 
Drunella  doddsi 
Cinygmula  sp. 
Epeorus  sp.-damaged 
Rhithrogena  sp. 


Cricotopus  (Isocladius)  Gr. 
Diamesa  sp. 
Micropsectra  sp. 
Orthocladius  sp. 
Tvetenia  sp. 


1 

0.32 

6 

CG 

21 

6.62 

0 

CG 

1 

0.32 

0 

CG 

17 

5.36 

4 

SC 

50 

15.77 

0 

SC 

14 

4.42 

0 

SC 

Total  Ephemeroptera 

104 

32.81 

Zapada  columbiana 
Megarcys  sp. 
Yoraperla  sp. 

46 

5 

14.51 
4.10 
1.58 

2 
2 
1 

SH 
PR 
SH 

Total  Plecoptera 

64 

20.19 

Parapsyche  elsis 
Rhyacophila  Hyalinata  Gr. 
Rhyacophila  narvae 

2 
'  1 

0.63 

10.73 

0.32 

1 

1 
1 

PR 
PR 
PR 

Total  Trichoptera 

37 

11.67 

Prosimulium  sp. 

1 

0.32 

3 

CF 

Total  Diptera 

1 

0.32 

3 

1.58 

93 

29.34 

1 

0.32 

4 

1.26 

3 

0.95 

7 

CG 

5 

CG 

7 

CG 

6 

CG 

5 

CG 

Total  Chironomidae 


106 


33.44 


Grand  Total 


317 


100.00 


Aquatic  Invertebrate  Summary  Data 


Site  Name:  Stillwater  River:  above  Daisy  Creek 

TOTAL  ABUNDANCE  ■       ,317 

Ephemeroptera  +  Plecoptera  + 

Trichoptera  (EPT)  abundance  205 

TOTAL  NUMBER  OF  TAXA  19 

Number  EPT  taxa  il2 

TAXONOMIC  GROUP  COMPOSITION 


Site  ID:  Y04USR01    8/2/01 


GROUP 

Misc.  Taxa 

Odonata 

Ephemeroptera 

Plecoptera 

Hemiptera 

Megaloptera 

Trichoptera 

Lepidoptera 

Coleoptera 

Diptera 

Chironomidae 


#TAXA       ABUNDAN  PERCENT 


1 

5 

1.58 

0 

0 

0.00 

6 

104 

32.81 

3 

64 

20.19 

0 

0 

0.00 

0 

0 

0.00 

3 

37 

11.67 

0 

0 

0.00 

0 

.  0 

0.00 

1 

1 

0.32 

5 

106 

33.44 

CONTRIBUTION  OF  DOMINANT  TAXA 

TAXON  ABUNDANCE 

Diamesa  sp. 

Epeorus  sp. -damaged 

Zapada  columbiana 

Rhyacophila  Hyalinata  Gr. 

Drunella  coloradensis 

SUBTOTAL  5  DOMINANTS 

Cinygmula  sp. 

Rhithrogena  sp. 

Megarcys  sp. 

Acari 

Yoraperla  sp. 

TOTAL  DOMINANTS 


SAPROBIC  INDICES 
HilsenhofTBiotic  Index 


rcE  1 

PERCENT 

93 

29.34 

50 

15.77 

46 

14.51 

34 

10.73 

21 

6.62 

244 

76.97 

17 

5.36 

14 

4.42 

13 

4.10 

5 

1.58 

5 

1.58 

298 

94.01 

2.55 


RATIOS  OF  TAX  GROUP  ABUNDANCES 
EPT/Chironomidae 


1.93 


FUNCTIONAL  FEEDING  GROUP  (FFG)  COMPOSITION 


GROUP 

Predator 

Parasite 

Collector-gatherer 

Collector-filterer 

Macrophyte-herbivore 

Piercer-herbivore 

Scraper 

Shredder 

Xylophage 

Omnivore 

Unknown 


RATIOS  OF  FFG  ABUNDANCES 
Scraper/Collector-filterer 
Scraper/(  Scraper  +  C.filterer) 
Shredder/Total  organisms 


#TAXA       ABUNDAN  PERCENT 


4 

50 

15.77 

1 

5 

1.58 

8 

129 

40.69 

1 

1 

0.32 

0 

0 

0.00 

0 

0 

0.00 

3 

81 

25.55 

2 

51 

16.09 

0 

0 

0.00 

0 

0 

0.00 

0 

0 

0.00 

81.00 
0.99 
0.0S 

DIVERSITY  MEASURES 
Shannon  H  (loge) 
Shannon  H  (log2) 
Evenness 
Simpson  D 


2.19 

3.17 
0.75 
0.15 


COMMUNITY  VOLTINISM  ANALYSIS 
TYPE  ABUNDANCE  PERCENT 

Multivoltine  85  26.89 

Univoltine  212  66.96 

Semivoltine  20  6.15 


#TAXA       ABUNDANCE  PERCENT 
Tolerant  1  1  0.32 

Intolerant  S  67  21.14 

dinger  12  164  51.74 


Aquatic  Invertebrate  Taxonomic  Data 


Site  Name:  Stillwater  River:  below  Daisy  Creek 


Site  ID:  Y04USR02  8/2/01 
Taxon 

Approx.  pcfxxnt  of  sample  used:  60 

Quantity          Percent 

HBI 

FFG 

Acari 

1                0.34 

5 

PA 

Total  Misc.  Taxa 

1                0.34 

Baetis  tricaudatus 
Cinygmula  sp. 
Epeorus  sp. -damaged 
Epeorus  grandis 


1 
3 

3 


0.34 
1.03 
1.03 
1.03 


CG 
SC 
SC 
SC 


Total  Ephemeroptera 


10 


3.45 


Plumiperla  sp. 
Sweltsa  sp. 
Zapada  columbiana 
Megarcys  sp. 
Yoraperla  sp. 


2 

1 

10 

1 

1 


0.69 
0.34 
3.45 
0.34 
0.34 


0 
1 
2 
2 
1 


PR 
PR 
SH 
PR 

SH 


Total  Plecoptera 


15 


5.17 


Parapsyche  elsis 
Micrasema  sp. 
Rhyacophila  grandis 
Rhyacophila  Hyalinata  Gr. 
Rhyacophila  Vagrita  Or. 
Rhyacophila  Vofixa  Gr. 


3 
2 

4 

10 


1.03 
0.69 
1.03 
6.21 
1.38 
3.45 


PR 
MH 
PR 
PR 
PR 
PR 


Total  Trichoptera 

40 

13.79 

Ceratopogoninae 
Oreogeton  sp. 
Rhabdomastix  sp. 

2 

.1 
1 

0.69 
0.34 
0.34 

6 
6 
3 

PR 
PR 
UN 

Total  Diptera 

4 

1.38 

Cricotopus  (Isocladius)  Gr. 
Diamesa  sp. 
Orthocladius  sp. 

181 

36 

3 

62.41 

12.41 

1.03 

7 
5 
6 

CG 
CG 
CG 

Total  Chironomidae 

220 

75.86 

Grand  Total 


290 


100.00 


Aquatic  Invertebrate  Summary  Data 


Site  Name:  Stillwater  River  near  wilderness  boundary 

177 
32 


Site  ID:  Y04USR04  8/24/01 


TOTAL  ABUNDANCE 
Ephemeroptera  +  Plecoptera  + 
Trichoptera  (EPT)  abundance 


TOTAL  NUMBER  OF  TAXA 
Number  EPT  taxa 

TAXONOMIC  GROUP  COMPOSITION 


14 

7 


GROUP 

Misc.  Taxa 

Odonata 

Ephemeroptera 

Plecoptera 

Hemiptera 

Megaloptera 

Trichoptera 

Lepidoptera 

Coleoptera 

Diptera 

Chironomidae 


#TAXA       ABUNDAN  PERCENT 


2 

79 

44.63 

0 

0 

0.00 

0 

0 

0.00 

2 

14 

7.91 

0 

0 

0.00 

0 

0 

0.00 

5 

18 

10.17 

0 

0 

0.00 

0 

.  0 

0.00 

1 

5 

2.82 

4 

61 

34.46 

CONTRIBUTION  OF  DOMINANT  TAXA 

TAXON  ABUNDANCE 

Acari 

Cricotopus  (Isocladius)  Gr. 

Polycelis  coronata 

Sweltsa  sp. 

Rhyacophila  Vofixa  Gr. 

SUBTOTAL  5  DOMINANTS 

Ceratopogoninae 

Megarcys  sp. 

Parapsyche  elsis 

Rhyacophila  Betteni  Gr. 

Rhyacophila  Hyalinata  Gr. 

TOTAL  DOMINANTS 


SAPROBIC  INDICES 
HilsenhofifBiotic  Index 


CE 

PERCENT 

61 

34.46 

55 

31.07 

18 

10.17 

10 

5.65 

7 

3.95 

151 

85.31 

5 

2.82 

4 

2.26 

4 

2.26 

3 

1.69 

3 

1.69 

170 

96.04 

4.79 


RATIOS  OF  TAX  GROUP  ABUNDANCES 
EIT/Chironomidae 


0.52 


FUNCTIONAL  FEEDING  GROUP  (FFG)  COMPOSITION 
GROUP  #TAXA       ABUNDAN  PERCENT 


Predator 

8 

37 

20.90 

Parasite 

1 

61 

34.46 

Collector-gatherer 

3 

79 

44.63 

Collector-mterer 

0 

0 

0.00 

Macrophyte-herbivore 

0 

0 

0.00 

Piercer-herbivore 

0 

0 

0.00 

Scraper 

0 

0 

0.00 

Shredder 

0 

0 

0.00 

Xylophage 

0 

0 

0.00 

Omnivore 

0 

0 

0.00 

Unknown 

0 

0 

0.00 

RATIOS  OF  FFG  ABUNDANCES 

Scraper/Collector-filterer 

#DIV/0! 

Scraper/(Scraper  +  C.filterer) 

#DIV/0! 

Shredder/Total  organisms 

0.00 

DIVERSITY  MEASURES 
Shannon  H  (loge) 
Shannon  H  (log2) 
Evenness 
Simpson  D 


1.84 
2.66 
0.70 
0.23 


COMMUNITY  VOLTINISM  ANALYSIS 
TYPE  ABUNDANCE  PERCENT 

Multivoltine  125  70.48 

Univoltine  41  23.31 

Semivoltine  11  6.21 


#TAXA  ABUNDANCE  PERCENT 

Tolerant  0                       0            0.00 

Intolerant  4                      16            9.04 

dinger  7                     77          43.50 


